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Abstract. It is known that, if a point in R n is driven by a bounded 
below potential V, whose gradient is always in a closed convex cone which 
contains no lines, then the velocity has a finite limit as time goes to +oo. 

The components of the asymptotic velocity, as functions of the initial 
data, are trivially constants of motion. We find sufficient conditions for these 
functions to be C k (2 < k < +oo) first integrals, independent and pairwise 
in involution. 

In this way we construct a large class of completely integrable systems. 
We can deal with very different asymptotic behaviours of the potential and 
we have persistence of the integrability under any small perturbation of the 
potential in an arbitrary compact set. 
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1. Introduction 



Given a smooth real function (p, q) h-> H (p, q) defined in an open domain fi of R n x R n 
we can consider the associated Hamiltonian system, that is, the autonomous system of 
ordinary differential equations 

. OH . OH 

q= ^p-> P = ~^q~- (L1) 

The function H is called a (time-independent) Hamiltonian. We remind that the Poisson 
brackets of two smooth real functions F, G: fi — >■ -R are 

{F ' G}:= S(£S~^©' (L2) 

An F G C 1 (0, R) is a constant of motion (or first integral) for the system (1.1) if and only 
if {F, H} = 0. 

Let us suppose that we find n functions Fi, . . . , F n : fi — > R of class C fc , 2 < k < +oo, 
such that: 

i) {F^ i7} = for all i (i.e., the Fi are first integrals of (1.1)); 

ii) {Fj, Fj} = for all i,j (i.e., the Fj are pairwise in involution); 

iii) VFi , . . . , VF n are linearly independent in all of O (the Fi themselves are then said to 
be independent). 

In this case a well known classical theorem says that the system (1.1) can be integrated 
by quadratures, in the usual sense of ordinary differential equations (see [2], Chapter 4, 
Section 1.1). 

If H itself is one of the functions Fi and the solutions of the Hamiltonian systems 
associated with each Fj are all global (i.e., defined on R), then the system (1.1) is called 
C k -completely integrable (see [2], Chapter 4, Section 1.2). Analogously we define analytic 
integrability. 

What is interesting about completely integrable systems is that the structure of the 
set of their solutions is very simple (see [2], Chapter 4, Section 1.2, Theorem 3). These 
properties, or rather the corresponding ones in a more general setting, are the foundations 
of a rich theory in the case when the level surfaces of the vector function (Fi, . . . , F n ) are 
compact. This paper does not deal with this last situation, but is concerned with proving 
the complete integrability of some systems with non-oscillatory behaviour, loosely related 
to scattering problems. 



Let V: R n — > R be a smooth function (called potential) and consider the Hamiltonian 
H(p, q) := ^\p\ 2 + V(q) with its associated system 

q = P, P = -VV(g). (1.3) 

Denote by t i— > (p(t,p,q),q(t,p,q)) the solution to (1.3) with (p,q) as initial data: 

p(o,P,q)=P, q(0,P,q) = q- 
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Our starting point is the following assumption on the potential V. 

Hypothesis 1.1 V is a function in C 2 (R n ;R) such that: 

i) V is bounded below; 

ii) there is a basis {61, ... , b n } for R n such that —W(q) ■ bi > for all q G R n and all bi. 
Of course, the system admits the first integral of energy 

l\p(t,P,q)\ 2 + V(q(t,p,q)) = 1 \p\ 2 + V(q) . (1.4) 

From i) we see that \p(-,p, q)\ must be bounded for each solution, so that by standard 
arguments in Ordinary Differential Equations we can prove that all solutions to (1.3) are 
defined for all times t G R. 

On the other hand, property ii) implies that t p(t,p, q) ■ bi is a monotone function 
for all (p, q) and for all b^. 

The whole of Hypothesis 1.1 thus ensures the existence, along each solution, of the 
following limit, the asymptotic velocity: 

Poo (p,q) := lim p(t,p,q)eR n . (1.5) 

t— >+oo 

The limit as t — > — 00 exists as well. 

These remarkably simple facts were pointed out by Gutkin in [5]. He called the 
potentials V satisfying Hypothesis 1.1 ii) cone potentials. The reason for this name is as 
follows. Let C be the convex cone in R n spanned by the forces — VV: 

C := I - A « VV(g a ) : ^ I finite set, A Q > 0, q a G R n Va G l\ . (1.6) 

and let V be the dual cone of C, defined by 

V := {w G R n : w ■ v > G C) . (1.7) 

Then Hypothesis 1.1 ii) means that T> has nonempty interior, or, equivalently, that the 
closure of C contains no straight lines (such cones C are called proper). We refer to Section 2 
for more details about cones. 

* * * 

Let us survey the content of the present paper. We are going to provide only hints to 
our assumptions and results. We will direct in each case to the precise statements scattered 
through the following Sections. 

In Section 2 we give a few generalities about cones in R n and prove a formula (Pro- 
position 2.4) that will be used extensively. 

Section 3 presents three simple instances of cone potentials for which the asymptotic 
velocity does not depend continuously on the initial data. The analysis of these counter- 
examples leads in Section 4 to write down our basic assumptions (the only global ones) 
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on the potential V, that, roughly speaking, amount to these: 

1) every level set of the potential V must be contained in a set of the form q + V, so that 
the asymptotic velocity turn out to belong to V (Hypothesis 4.1); 

2) the force — VV must push consistently toward the interior of T> ; somewhat less roughly, 
the component of — W(q) along any given direction of C shall be bounded below by 
a positive constant, when q varies on a (possibly noncompact) set of a certain sort 
(Hypothesis 4.2). 

Requirement 2) is actually the only severe limitation for our approach. In particular, it 
implies that C is contained in V, i.e., the scalar product of any two vectors from C is 
nonnegative (i.e., C has width not larger than 7r/2). Until Section 10 we will think of V as 
being defined on all of R n , but everything runs just as well if V is defined on a set of the 
form q + T>° . 

With the right hypotheses in hand, it becomes easy to prove that the asymptotic 
velocity always lies in the interior of the dual cone V (Proposition 4.3), with certain 
locally uniform estimates on the trajectories (Proposition 4.4). Such information is first 
used in Section 5 to find general sufficient conditions (Hypothesis 5.1) on the decay rate 
of V "at infinity" (in the direction of the cone V) for p x to be a continuous function of 
the initial data. The tools are the fact that p x can be expressed as an integral: 



and the theorems on uniform integrability. 

The first order differentiability of p x is less immediate. We get it in two different sets 
of assumptions. In Section 6 we impose an exponential decay on the second derivatives of V 
(Hypothesis 6.1). This will permit to exploit a simple Gronwall estimate on the solutions 
of the first variational equations of our system, and to use the theorems on differentiation 
under the integral sign in (1.8). In Section 7 we allow far more general asymptotics for V, 
but we add the side hypotheses of convexity on V and a kind of monotonicity in the Hessian 
matrix (Hypothesis 7.1). A Liapunov function built on the Hessian matrix of V will give 
a sharp control over the growth of the solutions of the first variational equation. As for 
the rest, Sections 6 and 7 run very much parallel to each other. Beside the mere regularity 
(Propositions 6.3 and 7.3), we also prove that p x , as a function of the initial data, is 
asymptotic, in the C 1 norm, to the projection (p,q) h-> p (Propositions 5.3, 6.4 and 7.4). 
This will be crucial in proving independence and involution in Section 9. 

In Section 8 we show how to get higher order differentiability of p x . This is not 
difficult, since the bulk of the job has already been done in Sections 6 and 7. 

In Section 9 we reap the rewards of the regularity theory to prove that the compo- 
nents of the asymptotic velocity are first integrals, independent and in involution, and to 
state the complete integrability of our systems (Theorem 9.1). Furthermore, we show that 
the potentials V satisfying our sufficient conditions for integrability can undergo arbitrary 
(small enough) perturbations on any compact set of R n without losing the property of 
yielding completely integrable systems (Persistence Theorem 9.2). The fact that the inte- 
grability is decided almost only on asymptotic behaviour and survives generic modifications 
in a bounded set seems to be unusual in the theory of integrable Hamiltonian systems. 




(1.8) 
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In Section 10 we give some examples. Namely, we provide manageable conditions 
(Hypotheses 10.1) on the functions /i, . . . , and on the vectors Vi, . . . , Vn (N > 1, no 
relation to n) so that our theory applies to the system with the potential V given by 

N 

V(q):=J2f<*(<l- v a)' 1 eRn (1-9) 

a=l 

(Proposition 10.5). A concrete instance is given in Corollary 10.6: if v a ■ vp > for all a, j3 
and if r > 0, then the Hamiltonian system with potential 

N 1 

V(g):=^ ? -, qe{qeR n : q-v a >0Va} (1.10) 

a=l 

is C°°-completely integrable. 

These cone potentials have polyhedrical (that is, finitely generated) cone C of the forces 
(Lemma 10.2). In a future paper (in preparation) we will provide an example where C is 
not polyhedrical. In fact, the present approach does not exploit such additional structures 
of V as being finite sum of one-dimensional functions. 

* * * 

An important analytically integrable system with cone potential (and cone wider 
than 7r/2) is the classical nonperiodic Toda Lattice system. It describes the dynamics 
of n particles on the line with coordinates qi, • • . ,q n interacting through an exponential 
potential 

n-1 

V(q) := » 4 = (91, • • • , Qn) G R n . (1.11) 

i=i 

This system was integrated by Henon and Flaschka (with different methods, unrelated to 
cone properties) — see the deep Lecture Notes [11]. 

Subsequently, Gutkin in [5] introduced cone potentials, recognizing that a large sub- 
class of the Toda-like potentials 

JV 

V(q) := c ^ Q ' Va , c Q >0, qeR n (1.12) 

a=l 

have the cone property (precisely, if the convex cone C spanned by the vectors v\ . . . , vn 
is proper). He noted that the components of the asymptotic velocity, which exists for all 
systems with bounded below cone potentials, are likely candidates for being n independent 
integrals in involution, as required by the definition of integrability. However this is not 
true in general, since they may even be discontinuous — see Section 3. 

In the case when C has amplitude not larger than tt/2, the conjecture (namely, C°°- 
integrability) was rigorously proved by Oliva and Castilla in [12] for the potentials (1.12) 
and also for (1.9), but only with functions f a having exponential asymptotic behaviour at 
+oo. Their method uses the finite-sum form of the potential to define a "compactifying" 
change of variables. Then they prove and apply a Lemma in Dynamical Systems (of inde- 
pendent interest too), concerning the differentiability of a foliation of invariant manifolds. 
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Oliva and Castilla drop the n/2 restriction for the following potentials: 



V(q) := e q - V0 + e' 



6 E R, qeR 3 , : 
i; fl = 0(O,l,-l) + (l 



0)v 



(1,-1,0), 



(1.13) 



with three degrees of freedom, and 




a > , q 



(q u q 2 ) G R 2 



(1.14) 



with two freedoms, and some generalizations thereof, all having either two or three free- 
doms. We think that admitting wide cones needs such strong restriction as low dimensions 
and/or specialized proofs exploiting the particular structure of a potential, and cannot be 
covered by a general theory of integrable systems with cone potential. The special role of 
the angle 7r/2 is not so surprising if one thinks about the behaviour of a billiard ball in a 
wedge. The dynamics is fundamentally simpler if the wedge is wider than tt/2, and this 
corresponds to a cone of the forces smaller than tt/2. 

As for the potential (1.14), Yoshida, Ramani, Grammaticos and Hietarinta [13], using 
Ziglin's methods [14], proved that the associated Hamiltonian system cannot be analyti- 
cally integrable if a 2 ^ m(m — l)/2 for m integer. Oliva and Castilla therefore remarked 
that these systems are C°° but not analytically integrable. For < a < 1 the system de- 
fined by (1.14) does fit into our framework too (Corollary 10.7). In particular, we cannot 
expect analytic integrability in the present approach either. 

Another well-known cone potential (also with cone wider than tt/2) yielding an ana- 
lytically completely integrable system is 



V ^) := ( . _ A2 ' qe{( qi ,...,q n )eR n : qi <q 2 <...<qn}. (1.15) 



It was introduced by Calogero and Marchioro (see [3], [8], and [4]) as the classical counter- 
part to a certain quantum mechanical system. Marchioro proved (among many other 
results) the integrability by explicit calculation in the case n = 3. The integrability in the 
general case was conjectured by Calogero and proved by Moser [10] through isospectral 
deformations. 

Moauro, Negrini and Oliva [private communication] introduced the potentials (1.10) 
with exponent r = 2 and put them into the framework of cone potentials. They proved the 
C°°-integrability for n = 2 and 3, even when the cone C has amplitude larger than tt/2. The 
compactification procedure, already successful in [12] for the exponential case (1.12), had 
to be supplemented with new ideas to overcome the degeneracies arising in this different 
situation. In particular, they use some interesting techniques that had been developed 
in [9] in connection with a Liapunov Stability problem. 

Acknowledgement. We thank Piero Negrini for a critical discussion on the manuscript. 



l<i<j<n 
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2. Cones 

Definition 2.1 A cone in R n is a nonempty subset C of R n such that v G C , A > 
Xv e C. 

All the cones we will consider are convex, and many of them closed, too. Note that 
the closure of a (convex) cone is also a (convex) cone. We denote the closure of a subset 
A of R n either as A or as c\(A). 

Definition 2.2 Given a convex cone C in R n , the dual of C is the set 

C* :={weR n : w ■ v > Vv G C} . 

It turns out that C* is a closed convex cone and that (C*)* = C. 

Definition 2.3 A convex cone C in R n is called proper if its dual C* has nonempty 
interior. 

It is easy to see that a convex cone is proper if and only if its closure contains no 
straight lines. In fact, C contains a straight line iff it contains both v and —v for some 
v 7^ 0. And this is equivalent to C* C {w : w ■ v = 0}. 

There is a neat way to express the distance of a point of a convex cone from the 
boundary in terms of the dual cone. We will repeatedly use this formula. 

Proposition 2.4 Let C be a convex cone in R n , not reduced to {0}, and let D be 
its dual. IfweD then 

dist(tu, dD) = min w ■ v . 

vec 

\v\ = l 

Proof. The distance of w form dD is the same as the distance from cl(R n \D) : since 
w G D. On the other hand, D is a closed convex cone, so we have 

c\{R n \D) = (J {ueR n : u ■ v < 0} . 

\v\=l 

In fact, from Hahn-Banach Theorem, for each u G c\(R n \D) : there is an affine function 
z i— > z ■ v + a, with \v\ = 1, that is nonnegative on the convex D and nonpositive on the 
half-line {6u : 9 > 0}. Since belongs to both sets, the constant a is zero. Since z ■ v > 
for all z G D, we have v G D* = (C*)* = C. This proves the inclusion C. On the other 
hand, if u G R n and there exists v G C, \v\ = 1, with u • v < 0, then u = lim-u n , where 
u n ■ v < 0. Hence u n G R n \D and «Gcl (R n \D). 
We can write 

dist(w,dD) = inf dist(to , {u : u ■ v < 0}) . 

i«i=i 

This last distance is the distance of w from a half-space. It thus coincides with w ■ v. The 
infimum is finally a minimum because the set {v G C : \v\ = 1} is compact and v w ■ v 
is continuous. 
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In the case of a polyhedral cone, the formula becomes a minimum over a finite set. 
We will use it in Section 10. 

Proposition 2.5 Let C be the closed cone generated 
R n \{0}. Then the dual cone D := C* is given by D = {w 
w G D we have 

dist(to, 3D) = minio • — ^— . 

a \v a \ 

Proof. The formula for D is easy. Therefore R n \D = \J a {u : u ■ v a < 0}. The 
union being finite, we can simply take the closure and write 

cl(R n \D) = [j{u : u ■ v a < 0} . 

a 

The rest of the proof is the same as for the Proposition 2.4. 



Remark 2.6 Let D be a convex cone in J? n ,with nonempty interior and not coincident 
with all of R n . From Hahn-Banach Theorem, we can separate the origin € R n from the 
interior of D, i.e., there exists v G R n \{0} such that v ■ w > Vmj G D°. 



by the vectors v±, . . . , vn G 
: w ■ v a > Va} and for all 
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3. Counterexamples to Continuity 

Consider the system (1.3) for the following potentials. 

3 

Counterexample 3.1 Let V: R — > R , q *-> e~ q . 

Obviously, this is a cone potential, that is, the conditions of Hypothesis 1.1 hold. In 
fact V > and C = R + , where C is the cone generated by the forces — see (1.6). Moreover 
there is the equilibrium (p, q) = (0, 0). 

For A > , let us consider (p(-, 0, A), q(-, 0, A)), i.e., the solution of (1.1) with p = and 
q = A as initial conditions. Then q(t, 0, A) — > +oo as t — > +oo for any A > . 
The first integral of energy (1.4) gives (for any t) 

\p(t,0,X)\ 2 + V{q(t,0,X)) =e~ x3 . 

Therefore there is a discontinuity of the asymptotic velocity (p,q) i-> p^p, <?) at (p,q) = 
(0, 0) since |Poc(0, A) | 2 — > 1 as A — > instead of = [p^O, 0)| 2 as for the equilibrium. 



In the previous counterexample the lack of continuity is related with the presence 
of an equilibrium. However, the absence of equilibria is not sufficient to guarantee the 
continuity, as the following example shows. 

Counterexample 3.2 Let V: R — >■ R , <?•->■— arctan q . 

It is a cone potential because V > — 1 and C = R+ . Furthermore, there are no 
equilibria. If the initial position is q = , and we conveniently choose the initial velocity 
p < 0, then the corresponding solution (p(t,p, 0), q(t,p, 0)) — > (0, — oo) as t — > +oo . From 
the conservation of energy 

^\p{t,p, 0) | 2 - arctan q(t,p,0) = ^|p| 2 , 

we see that the good choice for the initial velocity is p = —y/ir. 

Now, let us reduce the initial speed, so that the motion reverses its direction at a certain 
time. We easily see that the solution (p(-, — y/n + A, 0), (/(■, —\pn + A, 0)), for any A > 0, 
has the following asymptotic behaviour: 

q(t, -v / 7r + A,0) ->■ +oo, \p(t, -v/tt + A, 0)| 2 ^ 2n - 2\y/n + X 2 as t ->■ +oo . 

The further limit as A — > 0+ proves the discontinuity of (p, q) >->■ p^ (p, g) at the point 
(p, q) = (— s/tt, 0) (which are initial data of a solution with a vanishing asymptotic velocity, 
as we saw above). 



Of course the preceding counterexample is possible because V does not go to +oo as 
q — > — oo, i.e., because we do not have a "barrier" in the direction opposite to the force. 
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So far we have seen that we must avoid the equilibria and we need a "barrier" in order 
that any motion eventually "points in the sense of the forces". We shall give a precise 
formulation of these concepts in the next Section. Now let us give a last example to show 
that the barrier is not yet sufficient. For this we need two degrees of freedom. 

Counterexample 3.3 V: R 2 -> R , (q u q 2 ) H> e _<? i + e~ q2 . 

This is a cone potential, {(v 11 v 2 ) : v t > , v 2 > 0} U {(0, 0)} being the cone C of the 
forces, and the dual cone coinciding in this case too with C. We do not have equilibria. The 
behaviour of the solutions can be easily investigated because the two degrees of freedom 
are separate. By Counterexample 3.1, we see at once that a discontinuity in p x arises at 
the origin. 



What seems to go wrong in the third counterexample is that the force — W(q) does 
not drive toward the interior of the dual cone for the q along the axis q x = 0. 
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4. Geometrical Bounds for the Asymptotic Velocity 



The basic ingredient of this work is the potential function V, about which we started off 
with Hypothesis 1.1. From V we constructed the Hamiltonian function "H, the Hamiltonian 
system (1.3), its solutions (p(t,p, q), q(t,p, q)), the asymptotic velocity p^ip, <?), the convex 
cone C spanned by the forces (1.6) and its dual V := C* (1.7). 

The next assumptions on V are the first steps toward integrability. 

Hypothesis 4.1 For each M > there exists a q M E R n such that 

qeR n \(q M +V) V(q)>M. 

Hypothesis 4.2 For each q' , q" G R n such that q" G q' + V, and for each v G C\{0} 
there exists e > such that 

i^qEq' + V and q ■ v < q" ■ v ) =4> - VV(g-) • v > e . 

Hypothesis 4.2 implies in particular that — W(q) ■ v > for all q G R n , and all 
v G C\{0}. Therefore C C V and C has amplitude not larger than ir/2. 

Proposition 4.3 If Hypotheses 1.1, 4.1 and 4.2 hold, then 

p oo (p,q)-v>0 (4.1) 

for all v G C\{0} and all initial data (p, q). This is the same as saying that Poc(p, q) hes in 
the interior of the dual cone T>. 

Proof. Let us fix the initial data (p, q). The potential is bounded above along the 
trajectory: 

V(q(t,p,q))<^\p\ 2 +V(q) := M W G R (4.2) 

by the conservation of energy. Hypothesis 4.1 alone guarantees that q(t,p,q) remains in 
q M +V for all times. This already implies that p^ip, q) belongs to the closed set V. 
Let v G C\{0}. As we remarked immediately after Hypothesis 4.2, 

p(t,p,q) ■ v = -W(q(t,p, q))-v>0 

and so t i — y p(t, p,q)-v is an increasing funtion. 

We argue by contradiction. If p x (p, q) ■ v happened to be nonpositive, then 

<i(t, Pi Q) • v = p(t, p, q) ■ v < Vt G R 
and t i-> q(t,p, q) ■ v would be decreasing. Hence 

t>0 =^ q(t,p,q) • v < q • v . 
12 



Hypothesis 4.2 now yields that for t > the scalar product — VV(g(t,p, q)) ■ v is not less 
than some e > 0. Thus we can write: 

p(t,p,q)-v=p-v + / — W(q(s,p, q)) ■ v ds > 

Jo 

> p ■ v + et — > +oo as t — > +oo , 

and this contradicts the assumption p x {p, q) • v < 0. Finally, formula (4.1) is equivalent to 
Pcx>(Pi q) £ because of Proposition 2.4. 



In the sequel we will need the following information about the trajectories: locally 
uniformly in the initial data, 

1) the velocity p{t,p,q) is eventually in the interior of the dual cone, its distance from 
the boundary remains larger than a positive number 7, and 

2) the position q(t,p, q) enters and no longer quits any set of the form q + V. 

Proposition 4.4 In the hypotheses of Proposition 4.3, for each (p , q ) e R n x R n 
and each q e R n there esist 7 > 0, t E R and a bounded neighbourhood U of (p 0: q ) in 
R n x R n such that, for all t > t and (p, q) G U we have 

p(t,p,q)eV°, dist {p(t, p, q),dV)> 1 , 

q(t, p, q) Eq +V, dist (q(t, p, q),q + dV) >"f(t-t ). 



Proof. Since p^po^qo) G T>°, let 7 := (l/2)dist(p oo (p , ^ ), dV) > 0. Because 
of Proposition 2.4 and the continuity of the distance function, there exist t x G R and a 
bounded neighbourhood U of (p , q ) such that 

p(t 1 ,p,q)eV°, dtet(p(t lt p t q) t &D) > 7 

for all (p,q) G U. But for all (p, q) and all v G C\{0} the function t 1— >■ p(t,p,q) ■ v is 
increasing, so that the velocity p(t,p, q) lies in V° for all t > t 1 and all (p, q) G U, and its 
distance from the boundary is not less than 7. For all (p, q) G U, t >t 1 , v G C, \v\ = 1, we 
have 

(q(t, P, q) ~q ) -v = (q(t ,p,q) - q ) -v+ / p(s,p,q) ■ v ds > 

Jt\ 

> inf inf (q{t u p,q) - q ) ■ w + ^(t - := a + ^{t - U) , 

(p,q)eU -ec 

and finally, for t > t := max.{t 1 ,t 1 — (a/7)} the point q(t,p, q) belongs to q + V and 

dist( q(t,p,q),q + dV) = mm (q(t,p,q) - q ) -v > i(t-t ). 

M=i 
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We may ask what happens of V and VV along the trajectories. 

Proposition 4.5 In the hypotheses of Proposition 4.3, for all initial data (p, q) G 
R n x R n we have 

lim V(q(t,p, q)) — inf V . 

t— y+oo 

Proof. Fix e > and pick q e G R n such that V{q e ) < inf V + s. Let q G q s + V. 
Then 

V(q)-V(q £ ) = J ^V(q £ +e(q-q £ ))de = 

= [ VV(q s + 0(q-q £ ))-(q-q £ )d9<0, 
Jo 

because —W{q e + 9{q — q e )) G C and q — q £ G V = C*. Hence we can write 

qeq e + V V(q) <infV + £. 

On the other hand, Proposition 4.4 guarantees, in particular, that for all (p, q) there exists 
t £ G R such that 

t>t e q(t,p,q) Eq e +V. 

This concludes the proof. 

Corollary 4.6 In the hypotheses of Proposition 4.3, the following identity holds: 
n(p,q) = \\P\ 2 + V{q) = ^\ Poo (p,q)\ 2 + mfV. 

Within the assumptions of this Section, the gradient VV does not need to be infinites- 
imal along the trajectories. Already in one dimension, it is not difficult to figure out a 
V G C 2 (R, R) that decreases from +oo to 0, and whose graph has infinitely many smooth, 
but steep steps (whose height will obviously tend to zero): 

infV = 0, supV = +oo, V < , liminf V{q) < . 

q—y + oo 

All the q(t) will go to +oo as f -> +oo because our hypotheses are verified, so that they 
will never stop undergoing jerks (VV does not converge). The assumptions of the next 
Section will rule out this possibility. 
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5. Continuity 



This Section deals with the continuity of the asymptotic velocity with respect to the 
initial data. Gutkin in [7] already studied the problem, but in his setting he had no 
guarantee that p x belonged to the interior of the dual cone for all the trajectories. Much 
less did he obtain such crucial estimates as the ones in Proposition 4.4. So he obtained the 
continuity in a nonempty, open subset of the space of the initial data, defined in terms of 
p x itself. 

In our assumptions, we get global continuity. We will also prove an asymptotic prop- 
erty of p x , that will enable us later to determine the exact range of the mapping p x 
(Proposition 7.5). 

The asymptotic velocity p^ip, q) = lim^+oo p(t, p, q) can be expressed in terms of an 
integral: 



In the hypotheses of the last section, we know that q(-,p, q) is always contained in q M + T> 
(see formula (4.2)). Moreover, p^ip.q) is in the interior of the dual cone V, so that the 
distance of q(t,p, q) from the boundary of q M + V grows linearly as t — > +oo, as we saw 
in Proposition 4.4. We may expect p x to be a continuous function of (p, q) if the norm 
|VV(<z)| is dominated by an integrable function of the distance between q and q M + dV. 
We may thus use a uniform integrability theorem on the integral (5.1). 

Hypothesis 5.1 There exist q G R n and an weakly decreasing, integrable function 
h : R+ — > R such that 



Proposition 5.2 If Hypotheses 1.1, 4.1, 4.2 and 5.1 hold, then p x is a continuous 
function of the initial data. 

Proof. Let (p , q ) be fixed and pick 7 > 0, t E R and U from Proposition 4.4: 
q(t,p,q) eq +V, dist (q(t, p, q), q + dV) > j(t - t ) 
for all t > t , (p,q) € U. Using Hypothesis 5.1, 




(5.1) 





\VV{q{t,p,q))\ < h (dist(q{t,p,q),q + dV)^ < h ( 7 (t-t Q )) , 



so that we can apply the theorems on uniform integrability to the formula 




t 



and obtain our continuity result. 
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We are now provided with n continuous integrals of motion: the components of the 
asymptotic velocity p^ip, q). 

Roughly speaking, if we find a region q + T> where the driving force — VV is utterly 
negligible, we may expect that, if we start the motion there, with a velocity p in the interior 
of V, then the motion has approximately constant speed: 

p(t,p, q)^p 

so that, for those initial data (p, q) we have 

Poo(p,(?)~P- (5.3) 



Proposition 5.3 In the hypotheses of Proposition 5.2, for each /i > and each 
7 > 0, there exists q' G R n such that 

(pGX?°, dist(p,<9£>)> 7 , qeq' + V^j \ Poo (p, q) - p\ < // . 



Proof. Let \x > 0, 7 > be fixed and pick q from Hypothesis 5.1. Let d > be 
such that 

r+00 

/ h {^t) dt < ji . 

Choose q' £ q + V such that dist^, q + dV) > d . Then, for all p £ T>° such that 
dist(p, dV) > 7 and q E q' + V we have 



dist(q(t,p, q),q + dV) > -yt + d 



and so 



\p°°(p,q) -p\ 



-+00 



-WV(q(t,p, q))dt 



< 

r+00 r+00 
< / h (it + d )dt= / K('jt) dt < p . 

JO J d 
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6. First Order Differentiability 
without Convexity on the Potential 

In the Hypotheses 1.1 we know that the velocity p(t,p,q) is a differentiable function 
of the initial data at all finite times t. If we differentiate the equation 

P{t,p,q) =p(t ,P,q) + -W(q(s,p,q))ds (6.1) 

J t 

with respect to an arbitrary component of p or g, we obtain 

Dp(t,p,q) = Dp(t ,p,q) + -HV(q(s,p,q))Dq(s,p,q) ds , (6.2) 

Jto 

where we denote by HV the Hessian matrix of V and by D the partial derivative (we 
reserve the character D for the Jacobian matrix). 

In the hypotheses of Section 4, we know that formula (6.1) holds with t = +oo and 
the integrability is uniform. Can we expect the same for (6.2)? What we seem to need is: 

1) an a priori bound on the growth of Dq(t,p, q), locally uniform on (p, q); 

2) a rapid enough decrease of the norm of the Hessian HV along the trajectories t >->■ 
l(t,p, q). 

If the two estimates match appropriately, we can use the theorems on the differentiation 
under the integral sign. 

The function z(t) = Dq(t,p,q) satisfies the linear differential equation 

z(t) = -HV(q(t,p,q))z(t), (6.3) 

that can also be rewritten as a first-order system: 

lfz(t)\ = f In\fz(t)\ 

dt\z(t)J \-HV(q(t,p,q)) Oj\z(t)J- 

Let us denote by R(t, s,p, q) the evolution operator of the system, i.e., the 2n x 2n matrix 
solution of 

^R(t,s,p,q) = (-HV(q°(t,p,q)) 0* ) R ^ S,P > ^ ' R{s, s,p,q) = I 2n ■ 

Let n and n' be the two projections R n x R n — > R n defined as U(x, y) = x, U'(x, y) = y. 
Since 

Dq( u M) = mt , o, fi « (SSJJ;* J) = w*, o,a * (Z) ■ 

what we are interested in is the behaviour of ||ILR(£, s,p, q)\\ as t — > +oo. 
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We will carry out this program in two sets of hypotheses. In the rest of this Section 
our assumptions will be as follows. 

Hypothesis 6.1 There exist q 1 G R n , A 1 > 0, A x > such that 

qeq 1 +V \\HV(q)\\<A 1 exp(^-X 1 dist(q 1 q 1 + dV)Y 

On one hand, the mere fact that the Hessian is infinitesimal along the trajectories 
t i — v q(t,p, q) ensures, via a Gronwall lemma, that Dq(t,p, q) must grow less than exponen- 
tially as t — > +oo (i.e., it is o(e st ) for all e > 0). On the other hand the actual exponential 
decrease of the Hessian compensates for the other growth and yields the uniform integra- 
bility of (6.2). 

Lemma 6.2 Suppose that Hypotheses 1.1, 4.1, 4.2 and 6.1 hold. Then, for all e > 0, 
and for all t > 0, x, y G R n we have 

(peT>, q E q 1 + V , distfo q, + dV) > -y In 

=► |lli*(f , 0, p, q) Q I < 1 ((e|x| + \y\)e £t + (e\x\ - \y\)e- £t ) . 

Proof. Choose e > 0, q G q x + V such that dist(g, q 1 + dV) > (-1/AJ ln(e 2 /A 1 ). 
For any p G X>, we have 

dist (q(t,p, q),q! + dV) > dist(g, q 1 + dV) 

and so 

||HV(g(*,p,g))|| < A 1 e W (-\ 1 dist(q(t,p,q),q 1 +dV)^ < 

< A 1 exp^-A 1 dist(g,Q 1 + dvfj < e 2 . 

For any x, y G -R n , the function z(t) = UR(t, 0,p, q)(y) is a solution of (6.3), that can be 
rewritten in integral form: 

z(t) = x + ty + [ -(t- s)HV(q(s,p,q))z(s)ds. 
Jo 

Taking the norms, 

\z(t)\<\x\+t\y\+ [ (t- S )\\HV(q(s,p,q))\\\z( S )\ds< 
Jo 

<\x\ + t\y\ + e (t- s)\z(s)\ds. 
Jo 

A standard Gronwall argument yields that \z(t) \ < ip(t) for t > 0, where (p is the solution 
of 

(p(t) = \x\ + t\y\ + e / (t — s)(p(s) ds , 
and is precisely the expression appearing in the statement of the Lemma. 
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Proposition 6.3 Suppose that Hypotheses 1.1, 4.1, 4.2, and 6.1 hold. Then the 
asymptotic velocity p x is a C 1 function of the initial data. 

Proof. What we need is local uniform integrability of the integral in (6.2) for some 
t G R. Choose an initial condition (p ,q ). From Proposition 4.4, there exist 7 > 0, 
£i G R and a bounded neighbourhood U of (p , q ) such that 

p(t,p,q)eV°, dist(p(t,p,q),dV) > 7 , 

q(t,p,q) eq,+V, dist(q(t,p,q),q 1 + dV) > ^{t-Q 

for all t > ti, (p, q) G U. So we have 

\\W(q(t,p,q))\\ < Ae-^-^ 

for all t > £1, (p, <?) G Now choose £ > and £ > t x such that 

< e < A X 7 and A 1 e" Al7( * 0_tl) < e 2 . 

Since 

Dq(t,p, q) = UR(s, 0,p, q) =UR(s-t , 0,p(t o ,p, q),q(t ,p, q)) ^ 

(6.4) 

from Lemma 6.2 we get that 

\Dq(t,p,q)\ < ai e £{t - ta) + a 2 (6.5) 
for all t >t and all (p, q) G t/, where 

a x := sup — (e|.Dg(t ,p,(?)| + \Dp(t ,p,q)\) , a 2 := sup J|D^(t ,p, q)\ . 
(p,q)eu & v y (p,g)ec/ ^ 

We can finally write, for all t >t , (p, q) G U: 

I - W(q(t, p, q))Dq(t, p,q)\< A 1 e- Al7( *-* l) (a ie e(t - to) + a 2 ) 

and we are done. 



The approximate equality (5.3) extends to the derivatives of the functions involved. 
The character D stands for the Jacobian matrix. 

Proposition 6.4 In the hypotheses of Proposition 6.3, for each \i > and for each 
7 > there exists d > such that 



'peV°, dist(p,<9£>) > 7, 
q G <?! + V , dist((/, + <9£>) > d 



Dpoo(p,?) - Dn(p,g) 
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Proof. Fix 7 > 0. For any d > such that 

(Ae -A ldo) i/2 < Vy 

we set £ = (A 1 e _Ald °) 1 / 2 . Applying Lemma 6.2 we get that for all p E £>, q E q x + V such 
that 

1 e 2 

dist (g, q 1 + dV) > -—In — = d 

Ax A 1 

and for all t > 0, we have 

/nrr'X 1 1 

\Dq(t,p,q)\ = UR(t,0,P,q)l DU ) < -(s\DU'\ + \DU\)e st + -\DU'\ < 

< 1 ( ^1\DW\ + \DU\)e x ^ 2 + -\DU'\ . 
~ 2e V 2 1 / 2 1 1 

On the other hand, if moreover dist(p, dV) > 7, we have 

dist(g(t,p, q),q! + dV) > d + 7* , 

so that || HV(g(t,p, <?))|| < A 1 e _Al7t_Ald °. Putting the pieces together, and reminding that 

£ = (A.e-^) 1 / 2 : 



Dp 00 (p 1 q)-DU(p,q) 



+00 



-HV(q(t,p,q))Dq(t,p, q)dt 



< 



|-Dn|\ x 



2 ■ A l7 /" 2A l7 



It is clear that we can choose d so large that the last quantity is as small as we wish. 
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7. First Order Differentiability 
with convexity on the potential 

If we assume that the potential V is a convex function, then the Hessian matrix HV 
is nonnegative definite, and, from equation (6.3), 

z(t) = -HV(q(t,p,q))z(t), (7.1) 

it follows that z ■ z < 0. This lets us hope that we can derive a much sharper estimate on 
\Dq(t,p, q)\ than the mere less-than-exponential of Section 6. We will also assume that the 
quadratic form associated with HV((/) behaves monotonically with respect to q. Supposing 
V to be three times differentiate is not strictly necessary, but will simplify the proofs. 

Hypothesis 7.1 V is a C 3 function and there exists q 1 G R n such that 

i) V is convex on q 1 + V; 

ii) for all q' , q" G q x +T) and all z G R n we have 

q"eq' + V W(q")z ■ z < HV{q')z ■ z ; 

in) there exists a weakly decreasing function h x \ R + — > R such that J Q °°xh 1 (x)dx < +oo 
and 

qeq.+V \\HV(q)\\<h 1 (diat(q,q 1 + dV)y 

Lemma 7.2 Suppose that Hypotheses 1.1, 4.1, 4.2 and 7.1 i), ii) hold. Then, for all 
x, y G R n and t > we have 

/_ „ _ ^\ r |n^(t,o,^,^)(j)| < |^| + 1(!2/| + 1| hv(^)H 1/2 1^1) , 

r ' qe ^ +V ) {\U'R(t,0,p,q)Q\<\y\ + \\W^)\\ 1/2 \x\. 

Proof. Let p G V, q G q x + V be fixed. Let z(t) = UR(t, 0,p, q) Q), and consider the 
following Liapunov function: 

L(t, p, q) = \z(t)\ 2 + HV(q(t, p, q))z{t) ■ z(t) . 

L is a nonnegative quantity because V is convex in the points where the Hessian is evalu- 
ated. We are going to prove that L is decreasing in t for t > 0. Take the derivative with 
respect to t, and remind equation (7.1): 

j f L(t,p, q) = 2z(t) ■ z(t) + 2 W(q(t,p, q))z{t) ■ z(t) + (^( HV(g(f,p, q)) ))z(t) ■ z(t) = 

= \im J—(HV(q(s,p,q))- HV(q(t,p,q)))z(t)-z(t). 

This expression is nonpositive because of Hypothesis 7.1 ii) and because q(s,p,q) G 
Q(t,Pi q) + T> for s > t > 0. We thus have, for all t > 0: 

|i(t)| 2 < L(t,p,q) < L(0,p,q) < |i(0)| 2 + HV(q)z(0) ■ z(0) < \y\ 2 + HVfaJx - x , 

and so 

WR(t,0,p,q)(^j | = \z(t)\ < \y\ + || HV(^)I| 1/2 \x\ • 
The other inequality comes from the last one and from \z(t)\ < \x\ + J* \z(s)\ds. 





21 



Proposition 7.3 Suppose that Hypotheses 1.1, 4.1, 4.2 and 7.1 hold. Then the 
asymptotic velocity p x is a C 1 function of the initial data. 

Proof. Let (p , q ), 7 > 0, t t e R, U as in the proof of Proposition 6.3. Then, for 
all (p, q) e U and t > t t we have 

II HV( 9 (*,p,g))|| <M 7 (*-*0). 
On the other hand, from Lemma 7.2 and formula (6.4), we get that, again for all (p, 5) G £7, 

* > *i = 

where 

fll := sup \Dq{t u p,q)\, a 2 := sup (iDp^p, q)\ + || HV( ?1 )|| 1/2 \Dq{t^p, q)\) . 
(p,«)ei/ (p,g)€f/ v 7 

We can finally write, for all (p, q) E U, t > t x : 

\-HV(q(t,p, q))Dq(t,p, q)\ < (a, + a 2 (t - tj) h^t - tj) 



and we are done. 



Also Proposition 6.4 remains true in the modified hypotheses. 







Proposition 7.4 In the hypotheses of Proposition 7.3, for each fx > 0, and each 
7 > there exists d > such that 



peV°, dist(p,<9£>) > 7, 
qE q 1 +V , dist(g, ?! + <9Z>) > d 



Dp^p,^)- Dn(p,g) 



Proof. Fix // > and 7 > 0. Let p e V°, dist (p, dV) > 7, g e q 1 + V and 
dist (q, q 1 + &D) > d . Then, for all t > 

II HV(q(t,p,q))\\ <^(7* + do)- 
Moreover, from Lemma 7.2, we get 

|£>g(*,p,s)| = | ni? ( £ '°^^)(^n)l - l^n'l + ^pn'i + ll hv( ?1 )|| 1/2 |£>n|) := 

:= b 1 + b 2 t . 

In conclusion, 

r+00 

\D Poc (p,q)-DU(p,q)\= / - HV(q(t, p, q))Dq(t, p, q) 

Jo 



dt 



< 



< / (61 + b 2 t) h^t + do) dt . 
Jo 

It is clear that we can choose d so large that the last integral is as small as needed. 
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Proposition 7.5 In the hypotheses of Proposition 5.3 and either 6.4 or 7.6, the 
mapping p^ is surjective from R n x R n onto T>°. 

Proof. From Proposition 4.3, the image of p x is contained in T>°. To prove the 
reverse inclusion, let p G T>° . We can solve the equation p x (p, q) = p via the contraction 
principle. Let 7 := dist(p 0) dV) > 0. From Proposition 5.3 and either 6.4 or 7.6, there 
exists q G R n such that 

7 11 _ 11 1 

|p«»(p»?) -p\ < 2 ' II d Poo(p, «?) — Dn|| < - 

for all p G V such that dist(p, dV) > 7/2. Now the mapping p 1— > p + p — p x (ft, q) is a 
contraction of the closed ball {p G R n : \p — p \ < 7/2} into itself. The corresponding 
fixed point p solves p 00 (p, q) = p . Actually, we could make it without differentiability, in 
the mere hypotheses of Propositions 5.2 and 5.3, if we were willing to conjure up Brouwer's 
fixed point theorem. 
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8. Higher Order Differentiability 



Let us denote by D\, D 2 the partial derivative operators with respect to any two 
components of (p, q), by D 12 the second order derivative -D1-D2, and by D 3 V(g) the third 
differential of the potential V, regarded as a bilinear operator from R n x R n into R n , with 
the norm 

||D 3 V(<z)|| := sup \D 3 V(q)(x,y)\. 
N<i,li/I<i 

To get the second order differentiability of the asymptotic velocity p x , what we need 
is (local) uniform integrability of 



+00 







D 3 V(q(t,p, q)) (D iq (t,p, q),D 2 q(t,p, g)) - W(q(t,p, q))D^ 2 q(t,p, q)jdt. (8.1) 



To this purpose we must obtain an estimate on the growth of z(t) = D 12 q(t,p, q), which 
is a solution of the non-homogeneous linear differential equation: 

z(t) = -HV(q(t,p, q))z(t)-r(t,p,q), 
or, in first order system form: 

±(<t)\ ( 

dt\z(t)J \-W(q(t,P,d)) ; \z{t)) + \-r(t,P,q) 

Remind the evolution operator R introduced in Section 6. The function z(t) can be ex- 
pressed via R with the classical method of variation of the constants: 



(S)=/' ,i(M,ft »(-r(°ftJ' b <82> 



(note that in our case z(0) = i(0) = 0). 

We already have two sets of hypotheses that give an estimate of the evolution operator. 
All we are left to do is to give bounds on r(t,p, q). 

In the setting of Section 6, Diq(t,p, q) and D2q(t,p,q) grow less than exponentially 
as t — > +00. If we assume that || D 3 V(q(t,p, q))\\ decreases exponentially, then our scheme 
seems to work out. 

Hypothesis 8.1 V is a C 3 function and there exist q 2 G R n , A 2 > 0, A 2 > such 

that 

qeq 2 + V || D 3 V(q)\\ < A 2 exp(-\ 2 dist(q, q 2 + dV)^. 

We can safely assume that q 2 coincides with the q x of Hypothesis 6.1. 
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In the frame of Section 7, Diq(t,p, q) and D 2 q(t,p, q) grow linearly as t — > +00. 
Therefore the following assumption seems appropriate. 

Hypothesis 8.2 V is a C 3 function and there exist q 2 G R n and a weakly decreasing 
function h 2 : R+ — > R such that J + °° x 2 h 2 (x)dx < +00 and 

qeq 2 +V \\D S V(q)\\<h 2 (dist(q,q 2 + dV)y 



Here, too, we can assume q 2 to coincide with the (/! of Hypothesis 7.1. 

Proposition 8.3 Suppose that Hypotheses 1.1, 4.1, 4.2, 6.1 and 8.1 hold. Then the 
asymptotic velocity p^ is a C 2 function of the initial data. 

Proof. Let (p , q ), 7 > 0, U, e > 0, t x < t be as in the proof of Proposition 6.3. 
We showed there that 



\D iq (t,p,q) \ < a^-^ + a 2 , 



i = 1,2 



for all t > t and all (p, q) G U. This, together with Hypothesis 8.1 shows the uniform 
integrability of the first half of the integral (8.1): 

\r(t,p,q)\ < A^-^-^ { ai e e{t - to) + a 2 ) < a 3 e (£ - A27) * . 

We must now estimate z(t): 





*A = + f R(M ( - u s 

zit)) \D lt2 p(t ,p,q)J J to \-r(s,p,q) y 



We have 



R(t,s,p,q) = R(t-s,0,p(s,p, q),q(s,p, q)) 
From Lemma 6.2 we get a constant a 4 such that 



(8.3) 
(8.4) 



UR(t,s,p,q)^j \ < a 4 |y|e e ^- s) 



for all (p, q) G U, t < s < t. Hence 







UR(t,s,p,q)[ )ds 
t V r \ s iPi ( l) 



< a,a A / e e(t - s) e (£ - A27)s rfs < a,e et 



to 



Estimating \D 12 q{t ,p, q)\ on U by a constant, we can write \z(t)\ < a 6 e et . The last step 
is: 

I - HV(q(t, p, q))D lt2 q(t, p,q)\< A.e'^-^ a 6 e £t , 
for all (p, q) G U, t > £ , and the proof is complete. 
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Proposition 8.4 Suppose that Hypotheses 1.1, 4.1, 4.2, 7.1 and 8.2 hold. Then the 
asymptotic velocity p x is a C 2 function of the initial data. 

Proof. Let (p ,q ), U, 7, t x as in the proof of Proposition 6.3. In the proof of 
Proposition 7.3 we saw that 

\Diq(t,p, q)\ < a x + a 2 (t - t x ) , i = l,2 

for all (p, q) e U, t > t 1 . This, together with Hypothesis 8.2, gives 

\r{t,p,q)\ < (a x + a 2 (t - tj) h 2 (^(t - t t )) , 

and half of the job is done. Using again formulas (8.3) and (8.4) (with t 1 instead of t ), 
together with Lemma 7.2, we get 

/ U'R(t,s,p,q)( ,°__,V S ^ / \r(s,P,q)\ds< 
J tl \-r(s,p,q)J J tl 

r+00 

< / (a 1 + a 2 (s-t 1 ))h 2 ('y(s-t 1 ))ds :=a 3 , 
\z(t)\ <a 4 + a 3 := a 5 , \z(t)\ < a 6 + a 5 (t - t x ) . 

Finally: 

\-HV(q(t,p,q))D ia q(t,p,q)\ < (a 6 + a 5 (t - tj) h a (>y(t - tj) 
for all (p, q) e U, t > t x . The proof is complete. 



The hypotheses that guarantee higher order derivatives of p^ are now easy to guess. 
Denote by D m V the m-th differential of V, viewed as a multilinear operator form (R n ) m ~ 1 
into R n , endowed with the norm 

||D m V(g)|| :=sup{| D m V(Q)(x (1) ,...,x (m " 1) )| : x {i) G R n , \x {l) \ < 1 }. 
Hypothesis 8.5 

H m ) V is a C m+1 function and there exist X m > 0, A m > and q m E R n such that 
qeq m + V || D m+1 V(q)\\ < A rn exp(^-X rn dist(q,q rn + 8V)^j ; 

H' m ) V is a C m+1 function and there exist q m £ R n and a weakly decreasing function 
h m :R + — > R such that J ( ^ 00 x m h m (x)dx < +00 and 

qeq m +V \\D m+1 V(q)\\ < h m (dist(q,q m + 8V)Y 

The following proposition can be proved by induction on m, with essentially the same 
reasoning used in Propositions 8.3 and 8.4. 

Proposition 8.6 The asymptotic velocity p x is a C m function of the initial data, 
m > 2, if we assume Hypotheses 1.1, 4.1, 4.2, and either i) or ii) of the following: 

i) Hypothesis 6.1 plus H 2 , H 3 , ... , H m of Hypothesis 8.5; 

ii) Hypothesis 7.1 plus H' 2 , H' 3 , ... , H' m of Hypothesis 8.5. 
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9. Complete Integrability and Persistence 



It is time to exploit the regularity theory developed so far to achieve our main goal: 
the integrability of the system 

?=7^> ^ = "^' nP,Q)--=^\p\ 2 + y(Q)- (9-1) 

We introduce the notation Xf to mean the Hamiltonian vector field defined by some 
smooth /: R 2n — > R. In particular is the vector field in (9.1). We will say that Xf is 
complete if all the solutions of the Hamilton equations q = df /dp, p = —df/dq are global, 
that is, defined on the whole of R. 

We are going to prove that the components of the asymptotic velocity are n first 
integrals independent and in involution. Moreover we can include the Hamiltonian H into 
a set of n globally independent first integrals in involution T\, . . . , J- n , whose associated 
vector fields Xj x , . . . , Xj n are complete. The J^, 2 < % < n, will be obtained from p x 
through a linear transformation. We will use the fact that, in our hypotheses, tjIPoJ 2 is 
just "H, up to an immaterial additive constant (Corollary 4.6). 

Theorem 9.1 (Complete Integrability) Assume the hypotheses of Proposition 
either 8.3 or 8.4. Then the n components 

Poo,l 5 • • • 5 Poo,™ 

of the asymptotic velocity are independent C 2 integrals of motion and they are (pairwise) 
in involution. This means that, for all (p, q) G R n x R n , the gradients 

VlVxfe q) , VPooAPi VPoo, n (p, q) (9-2) 

are linearly independent and the Poisson brackets vanish identically: 

{Poo,i»Poo,j}(P,9) = 0. 

Hence the system (9.1) is integrable by quadratures. 

Furthermore, there exists an orthogonal transformation A: R n — > R n such that the 
functions T\ , . . . , T n defined as 

^ " H (9.3) 
Ti := {Apjji for i = 2, 3,..., n, 

are independent C 2 integrals of motion, (pairwise) in involution and the Hamiltonian vector 

fields Xj7. (1 < i < n) are complete. 

In particular, the system (9. 1 ) is completely integrable. 
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Proof. Claiming that the gradients in (9.2) are independent is the same as saying 
that the Jacobian matrix Dp^ (p^ thought of as column vector) has maximum rank. 
Denote by {<& t } t eR the flow of the system (1.1), i.e., <fr*(p, q) = (p(t,p,q),q(t,p,q)). Since 
p^ is an integral of motion, 

Poo o VteR. 
By differentiating at (p, q) we have 

Dp^^ip, q)) D$*(p, q) = D Poo (p, q) . 

From this we have that 

rank Dp^ 

is a constant of motion, because D<£>*(p, q) is invertible. 

From Poisson's Theorem (see e.g. [1] , Section 40) we know that the Poisson brackets 

{Poo.i , Poo, j } 

are also constants of motion. 

We now use Proposition either 6.5 or 7.6. Along any trajectory, the velocity p{t,p,q) 
eventually enters T>° and keeps from its boundary a distance larger than 7 > 0. Moreover, 
q(t,p, q) enters all the sets of the form q +T> (Proposition 4.4). Hence, along any trajectory 
the derivatives of p x tend to the derivatives of the projection IT and we can compute: 

rank Dp^p, q) = rank Dp^f&faq)) = lim rank Dp x ($*(p, q)) = 
= rank DIT = n 

2 2 

(the set of the n x 2n matrices with maximum rank is open in R n ) and 

{Poo,,, ,Poo, 3 }(p,(Z) = {Poo.^ , Poo,,}($*(p,?)) = lim {p^, p 00 , j }($ t (p,g)) = 

= {n i ,n j } = o 

Now, let us prove the second part of the theorem. We can use the Remark 2.6 to 
construct an orthogonal transformation A : R n — > R n (we will write A also for the 
associated matrix) such that the first component (Aw)± is strictly positive for all w G V°. 
Hence, from Proposition 4.3 we have that 

(Ap oo (p,g)) 1 >0 V(p,g). (9.4) 

Define Ti, T2, • • • , T n as in (9.3). Proposition 4.5 implies that J-\=% = ^Ip^p+inf V. Let 
DJ 7 be the Jacobian matrix of T := (J-i, . . . , J r n ) T (the T means transposition of matrices). 
Proving that T\, . . . , T n are independent is equivalent to proving that the rank of DJ 7 is n. 
Since J 7 is a first integral of motion, the rank of DJ 7 does not change along each trajectory. 
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From Proposition either 6.5 or 7.6 we know that along all trajectories the matrix Dp^ 
tends to the matrix DII = (J n ,0), so that D(Ap 00 ) ->■ D(AU) = (A, 0) = (/ n ,0)(g), and 
DJ 7 ! tends to (p^,0) = ((^Ap^) 7 , 0) = ((A Poo ) T , 0) (£) . Thus we can write 



rank DJ 7 (p,q)= lim rank DJ 7 (^ t (p,q)) 

t— >+oo 



rank 



/ (^Poo)l 

1 
1 



V 











{A Poa ) 7 



1 



n 



because of (9.4). So T\, . . . , T n are independent. 
We have 

n 

{u , j-j = {% , (A Poo ) j = ^{n , Poo ,,} = o 

3=1 

because the Poisson brackets are bilinear. Therefore, T\, . . . , T n are pairwise in involution 
since p Xil , . . . t p Xin are. 

The last result in particular says that T-L is a first integral of each Xjr. . This implies the 
completeness of each vector field X^. as we can see by the same argument which gave the 
completeness of X n in Section 1. 





We remark that we could have included the Hamiltonian T-L into a set of n locally inde- 
pendent integrals of motion, extracted from p x , in the following way. We could have simply 
chosen a nonvanishing component of p x (there is always one, locally, because p^ G V°) 
and replaced it with T-L. This would give independence generally only in a neighbourhood 
of each level set of p x , unless we have been lucky or deft from the very beginning in the 
choice of the ortho normal reference basis of R n . 

It is legitimate to ask what happens to the integrability of our system if we perturb 
the potential V in a compact set K: 

H{p,q):=\\p\ 2 + V{q) + f{q), (9.5) 

with / a smooth function vanishing outside K. The new potential V + / does not need to 
be a cone potential, or, even if it is, it may not have the same cone. Think for example to 
the case when the cone C of V has empty interior. The global hypotheses 4.1 and 4.2 are 
very sensitive to the cones C and T>, and there is no hope to verify them for the new cone. 
Nevertheless, we will prove that if |V/| is sufficiently small, then all the trajectories of the 
new system eventually quit K for good, the potential V leads henceforth undisturbed, and 
all our conclusions about regularity and integrability remain true. 
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Theorem 9.2 (Persistence) Suppose that V verifies the hypotheses of Theorem 9.1. 
Let K C R n be compact. Then there exists an e > with the following property. If 
f: R n — > R is a C 3 function with support in K and 

sup|V/|<£, (9.6) 

then the Hamiltonian system whose Hamiltonian H is given by Equation (9.5) is completely 
integrable. Namely, all trajectories have asymptotic velocities for which all the claims in 
Theorem 9.1 apply. 

Proof. We still denote by C and V the cones associated with V. Let v G C\{0}. 
Since K is compact and V has nonempty interior, there exist q\ q" G R n , q" G q' + V, 
such that 

K C q + V , q-v<q"-v Vg G K . 

Define e by 

e := inf {— W(q) ■ v : q G q + V , q ■ v < q" ■ v} . 

From Hypothesis 4.2, e is positive. Suppose that / is a C 3 function with support in K 
and verifies (9.6). Denote by (p(t,p, q), q(t,p, q)) the trajectories corresponding to the new 
Hamiltonian (9.5). We claim that 

V(p, q) 3t G R such that Vt > t q(t, p,q) -v > q" -v . (9.7) 

In fact, there certainly exists q G R n such that q(t,p,q) G q + V for all t G R because 
Hypothesis 4.1 still holds for V + / and V. We can safely assume q' G q + V. Let 

e' := inf{-VV(g) • v : qEq + V, q ■ v < q" ■ v) . 

We can write 

(qeq + V and q ■ v < q" ■ v ) -V(V + /)(?)• v > min{e - inf |V/| , e'} > . 

Now the same reasoning as in the proof of Proposition 4.3 leads to (9.7). 
To conclude, we only need to note that the t in (9.7) can be chosen locally independent 
of (p, g), and that in all the hypotheses on V after Section 4 we can always suppose q , 
q x ... to belong to q" + V. 
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10. Examples 



We are going to present a class of examples for which our assumptions for complete 
C 2 integrability hold, in the convex case. Conditions for C k integrability, 2 < k < +00, 
are easily added. 

Hypotheses 10.1 For a = 1, . . . , N, let f a be a C s real function, defined either on 
all of R or on the interval ]0, +00 [, such that, for all a, 

inf/ a = 0, sup/ Q = +oo. (10.1) 

f a (x)<0 Vx, (10.2) 

%(x)>0 Vx>l, (10.3) 

f£ is weakly decreasing on [1, +oo[ , f£{x) < Vx > 1 . (10.4) 
Let v u ...,v N G R n \{0} be such that 

v a -v/3>0 Va,/3. (10.5) 
Define the following potential V: 

JV 

V(?):=$] (10.6) 



a=l 



either on J? n or in the set {w G -R n : w ■ v a > Va = 1, . . . , iV} . The gradient, Hessian 
and third differential of V are given by 



JV 



VV(q) = J2f'a(<l-Va)v a , (10.7) 

a=l 
JV 

HV( 9 ) = ^ • V«) ® v a (10.8) 

a=l 

JV 

D 3 V(q) = ^ f™(q ■ v a ) v a ®v a ®v a , (10.9) 



a=l 



where ® indicates the tensor product: (u a ® fa)i,j : = Va,iV a ,j, etc. 

From (10.1), (10.2), (10.5) and (10.7) it is clear that Hypothesis 1.1 is satisfied for V. 
Also the other requirements for integrability are met, and the proof will be made in several 
steps, culminating in the statement of Proposition 10.3 and its corollaries. 
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It will be convenient to have from the start ag 6 R n such that 

q ■ v a > 1 Va (10.10) 

(for example, q := pY^p=i v P-> f° r P large). Recall Proposition 2.5 and note that, for all 
q e q + V, a = 1, . . . ,N 

q ■ v ao = {q~ q ) ■ v ao + q ■ v ao > min \vp\ dist(g, q + dV) + 1 . (10.11) 



Lemma 10.2 The closure of the convex cone C generated by the forces — VV coincides 
with the convex cone generated by the v a : 

N 



and the dua] VofC is given by {w : w ■ v a > Vet}. 

Proof. Denote by C the cone generated by v±, . . . , vn and by D its dual. Of course 
C C C, so that D d T>. Take to G <9-D. We are done if we show that w G <9£>. Consider 
the following two sets of indices: I\ := {a : w ■ v a > 0}, I2 := {/3 : w ■ vp = 0}. From 
Proposition 2.5 and the fact that w G dD we deduce that I2 7^ 0. Evaluate now — VV 
along the line q + tw (q given by (10.10)), for r — > +00; the terms f' a ((q + rw) ■ v a ) 
vanish for a G ii, so that 

C 9 -VV(g„ + ™) -> - ^ /£( 9o -i^)^ := u G (J. 

We have 6 7^ because —f' a > and the vp lie in the interior of a half-space (the cone C 
is proper) . From the definition of I2 we see that w ■ v = and finally w G <9£> because of 
Proposition 2.4. The formula for V is an easy consequence. 



Lemma 10.3 Let g: [0, +00 [— > R be a nonnegative, weakly decreasing function such 

that 

r+00 

/ x m g(x) dx < +00 
for some integer m > 0. Then lim^-^+oo x m+1 (7(x) = 0. 
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Proof. Suppose there exists Xi /* +00 such that x^ l+1 g(x i ) > e > for all i > 1. 
Then 

x m 

Xi-i < x < Xi x m g{x) > x m g(xi) > e—^ . 



X i 



Define h: [0, +00 [— > R as 

h(x) 



if0<x<xi, 
ex m /x 7 [ l+1 if Xi-i < x < Xi, i > 2. 



Then h(x) < x m g(x) and 

r+00 +00 m+l _ m+l +00, 

/ h(x) dx = V ■ - i ~ 1 = £ V 1- M ) 

y fe^r +1 m+l m+1 fev ^r +1 / 



We can safely assume, for instance, that Xi > 2xi-i for all i > 2, and this yields /i, and 
hence x m g(x) too, to have infinite integral, against the hypothesis. 



Lemma 10.4 Let f G C M ([0, +oo[), M > 1, and suppose that, for all x > 0, 
m = 0,...,M, 

f(m)(\(>0 ifm is even, 
1 lXJ l<0 if mis odd. 

Then / +oo x m - 1 |/( m )(x)|da;< +00 for all m = 1, . . . , M . 

Proof. It is certainly true if M = 1. Suppose it is true for M — 1 and write 

/ x M - 1 f( M \x)dx = x M - 1 f( M - 1 \x)\+™-(M-l) / x M - 2 / (M_1) (^)^. 
</o Jo 

The last integral converges for the induction hypothesis. The term x M ~ x f^ M ~ 1 \x) is 
infinitesimal as x — > +00 again because x M ~ 2 f^ M ~ l \x) is integrable, with the help of 
Lemma 10.3. 



Verification of 4.1 Let M > 0. From 10.1 and 10.2 we get x M G R such that 
fa(x)>M \/x<x M Va = l,...,iV. 

Let q M G -R n such that 

q M ■ v a < x M Va 

(for example, q M = 9J2 a v a for 6 negatively large). Then, for all q G R n \(q M + V) there 
exists a M such that q M ■ v aM < x M and hence 

V(<Z) > /a M (g • Va M ) > U M (X M ) >M. 
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Verification of 4.2 Let q' belong to the domain of V and q" G q' + V. Let v G C\{0}. 
Up to a reordering of indices, we can write 

N' 

v = Y^ C P V P e C with 1 < N ' < N and C I3 > • (10.12) 

0=1 

By (10.7), (10.12), (10.2) and (10.5) we have 

N N N' 

-VV(q) -v = -^2f' a (q-v a )v a -v = ^ ^ cp (-f' a (q • v a )) v a -vp> 

a=l a=l /3=1 

AT' 

0=1 

Moreover 

q e q' + V => V/3 g' • u p < q ■ vp . 

On the other hand 

q ■ v < q" • v =>- 3/3 () such that q ■ vp < q" • vp . 

Hence 

^ q G q' + V and q • v < q" • v ^ =>• 3/3 such that (/' • vp < q ■ vp < q" ■ vp . 

For each (3 = 1,..., N', define 

e p := mm{-cpfp(x) : q' ■ vp < x < q" ■ vp} > . 
We can conclude that 

(^qEq' + V and q ■ v < q" ■ v ) =4> -VV(<?) • f > mine^ > . 



Verification of 5.1 Recall q from formula (10.9). Since x \-¥ \f' a (x)\ is weakly 
decreasing on [1, +oo[ and from (10.11), we can compute, for all q G q + V, 

N 

|VV(g)| < J] ' 1^1 ^ ^o(dist(g,g + 9P) ) , 

a=l 

where 

AT 

/i (x) :=XKI |/4(nun|u j9 |x+ 1)| , 

a=l 

and /i is weakly decreasing and integrable on [0, +oo[ because each |/^| is weakly decreasing 
and integrable on [1, +oo[. 
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Verification of 7.1 i). From Formula (10.3) we see that each f a is convex on [1, +oo[. 
So the potential V is convex on q + V, q given by (10.10), because it is sum of convex 
functions. 



Verification of 7.1 ii). From Formula (10.4) we see that /" * s a weakly decreasing 
function on [l,+oo[. If q',q" G q +V, q" G q'+T>, then we have q" • v a > q'-v a > q -v a > 1 
and 

Hence, from (10.8), 

JV 

HV{q")z.z=Y J f'M'-v a ){v a -zf< 

a=l 
JV 

<Y,f'^ - v o){v a -z) 2 = HV{q')z-z. 

a=l 

Verification of 7.1 iii). We have, from Formula (10.8) and (10.11), for all q G q + V, 

N 

II HV(g)|| < ^2\fZ(q ■ v a )\ \v a \ 2 < h, (dist(g, q + dV) ) , 

a=l 

where 

JV 

^(x) := ^ |v c J 2 |/^(min|v /3 |:r + 1)| , 

a=l 

and h 1 is weakly decreasing on [0, +oo[ and J^°°x\h 1 (x) \ dx < +oo because of Lemma 10.4. 



Verification of 8.2 From (10.9) and (10.4) we have, for all q G q + V, 

II D 3 V((z)|| <fc 2 (dist(g,g + d£>)), 

where 

JV 

/i 2 (x) := \ v ct?\fa\™™\v a \x+ 1)| 
a=l 

and, as usual, /i 2 is weakly decreasing on [0, +oo[ and § Q °°x 2 h 2 (x) dx < +oo. 
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Proposition 10.5 Suppose that the functions /i, . . . , /jv & n d the vectors vi,...,vn 
satisfy the Hypotheses 10.1. Let the potential V be defined by Formula (10.6). Then 
Theorem 9.2 applies, so that the Hamiltonian system 

q=p, p=-W(g) 

is C 2 -completely integrable. 

If we assume, moreover, that the f a are C k+1 , 2 < k < +oo, and that for all 4 < m < k + l, 
a = 1, . . . , N, x > 1 we have 

f(m)/ \ j > if m is even 
J a K x ) |< o ifm is oddj 

(and that be weakly decreasing if k < +oo), then is C k . 

Corollary 10.6 Let vi, . . . , vn £ R n \{0} be such that v a ■ v p > for all a, /3. Let 
r > and define the potential 

on the set T>° = {q G R n : q ■ v a > Va}. Then the associated Hamiltonian system is 
C '°° -completely integrable. 

Corollary 10.7 Let v±, . . . ,vn £ R n \{0} be such that v a ■ vp > for all a, f3, and 
let c a > 0. Define the potential 

N 

V(q) :=^c a e-^ 

a=l 

on R n . Then the associated Hamiltonian system is C°° -completely integrable. 
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